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It is widely recognized that when classical optimal strategies are applied with parameters esti- 
mated from data, the resulting portfolio weights are remarkably volatile and unstable over time. 
The predominant explanation for this is the difficulty of estimating expected returns accurately. 
In this paper, we modify the n stock Black-Scholes model by introducing a new parametrization 
of the drift rates. We solve Markowitz' continuous time portfolio problem in this framework. 
The optimal portfolio weights correspond to keeping 1/n of the wealth invested in stocks in 
each of the n Brownian motions. The strategy is applied out-of-sample to a large data set. The 
portfolio weights are stable over time and obtain a significantly higher Sharpe ratio than the 
classical 1/n strategy. 
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1. Introduction 

The fundamental question of portfolio optimization - How do we trade in the stock 
market in the best possible way? - is as challenging today as ever. Classical strategies, 
such as Markowitz' mean-variance portfolio, applied with parameters estimated from 
data are known to give exceptionally volatile portfolio weights. This is primarily due to 
the difficulty of estimating expected returns with sufficient accuracy; see the examples 
in [2] . In this paper we develop a new approach to circumventing this difficulty. 

Several different methods for resolving this difficulty have been published. For example, 
Black and Littcrman [2] proposed to estimate the expected returns by combining Capital 
Asset Pricing Model (CAPM) equilibrium with subjective investor views. A drawback 
with this approach is that the investor's beliefs must be quantified by specifying numbers 
for both the expected returns and the uncertainty in them. 

The Arbitrage Pricing Theory (APT), see [11], is another acclaimed approach. The 
APT models the discrete time returns of the stocks as a linear combination of indepen- 
dent factors. The APT relies on statistical estimates of the expected returns that are 
constructed to fit historical data and hence again may lead to unstable portfolio weights. 
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Yet another popular method for dealing with the difficulty of estimating expected 
returns is simply to ignore them. This idea is pursued in the classical 1/n strategy which 
puts 1/n of the investor's capital in each of n available assets. Intuitively this strategy 
should be well diversified. However, this may not be the case due to covariation between 
different stocks. Since it is possible to obtain good estimates of covariances between stock 
returns, we want to use this information in our portfolio construction. 

Recently, it has been proposed to let expected returns depend on ranks. These ranks 
could, for example, be based on the capital distribution of the market, assigning rank 1 
to the stock with the highest market capitalization, rank 2 to the second highest, and so 
on. For developments of this idea, see [3] . 

Modern portfolio optimization was initialized by Markowitz in [8]. Markowitz mea- 
sured the risk of a portfolio by the variance of its return. He then formulated a one-period 
quadratic program where he minimized a portfolio's variance subject to the constraint 
that the expected return should be greater than some constant. Merton ([9] and [10]) was 
the first to consider continuous time portfolio optimization. He used dynamic program- 
ming and stochastic control to maximize the expected utility of the investor's terminal 
wealth. The first results on continuous time versions of the Markowitz problem were 
published rather recently; see [1, 5, 6, 7, 12, 13]. 

The goal of this paper is to find trading strategies that circumvent the problems associ- 
ated with estimating expected returns. Further, we want to include options for investors 
to specify their unique market views in a flexible and non-numerical way. To this end, the 
geometric Brownian motion is used as our n stock market model with the modification 
that the Brownian motions all have equal positive drifts. The classical assumption is to 
assume zero drift for all Brownian motions. Hence, we obtain an explicit connection be- 
tween risk and return by model construction, as the volatility matrix determines both the 
covariance matrix and the expected returns for the stocks. In effect, the expected returns 
for the stocks are determined by how each stock is exposed to the underlying Brownian 
motions. This connection is missing in the classical paramctrization of the Black-Scholcs 
model. Further, given a covariance matrix, the volatility matrix is non-unique. This al- 
lows the investor to impose views on the market by selecting an appropriate volatility 
matrix, which implies expected returns of the stocks that are consistent with those views. 

We solve Markowitz' continuous time portfolio problem explicitly for our n stock mar- 
ket model. The optimal strategy tt* corresponds to holding 1/n of the wealth invested 
in stocks in each of the n underlying Brownian motions. This is not the same as holding 
1/n of the wealth in each stock. 

We apply tt* , out-of-samplc, to a large set of daily market index data. The long- 
term performance of tt* is investigated under the assumption that the investor has no 
preferences or market views. We find that tt* is stable over time and outperforms all but 
one of the underlying market assets in terms of Sharpe ratios. Moreover, we can reject the 
hypothesis that the classical 1/n strategy gives a higher Sharpe ratio than tt* very clearly. 

We present our model in Section 2. Further, we indicate some procedures for obtaining 
volatility matrices that are consistent with investor specified market views. A continuous- 
time version of Markowitz' problem for n stocks is solved explicitly in Section 3. Section 
4 contains an empirical study of our strategy. 
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In this section we present the model for the stocks. Further, we discuss how to estimate 
the volatility matrix such that the expected returns for the stocks reflect the investor's 
market views. 



2.1. The stock price model 

For < t < T < oo, we assume as given a complete probability space (CI, J 7 , P) with a 
filtration {J-t}o<t<T satisfying the usual conditions. We take n independent Brownian 
motions Bi and define the stocks Si, i = l,...,n, to have the dynamics 

dSi(t) = Si(t) (r dt + J2 °i,3 [0* ~r)dt + dBj (t)]j , (1) 

where r > is the continuously compounded interest rate, the constant /x > r is a drift 
parameter, and the volatility matrix a := {&i,j}fj=i is assumed to be non-singular. The 
stock price processes then are 



-Y^^A(^-r)t + B 3 {t)]\, (2) 
i=i / 



for -Bi(O) = • • • = B n (0) = 0. We also equip the market with a risk-free bond with dynam- 
ics 

dR(t) = rR(t) dt. 

This is a parametrization of the geometric Brownian motion, known in finance as the 
Black-Scholes model, with the modification that the Brownian motions are assumed 
to have equal positive drifts fi — r. Note that the Brownian motions are not intended 
for direct interpretations as observable quantities. They are merely a mathematically 
necessary partition of the randomness to avoid arbitrage possibilities. 

We denote the covariance matrix by C and the continuously compounded rates of 
return for the stocks Si by i = 1, . . . ,n. Equation (2) gives that 

n 

ti = r + (v- r )^2^j, (3) 

3=1 



from which we can calculate the expected returns. Further, for each stock i, the market 
price of risk 
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so we can write 



Note that since the row sums of the volatility matrix determine the expected returns 
for the stocks are determined by how each stock is exposed to the underlying Brownian 
motions. We will sec later that the optimal strategy for Markowitz' problem in continuous 
time for this model can be applied without knowing the parameters fi and r, as long as we 
have the volatility matrix a. Hence we have partly circumvented the problems associated 
with estimating expected returns for the purpose of portfolio optimization. 
We next discuss how to choose the volatility matrix. 

2.2. The volatility matrix and the rates of return 

In mathematical finance, the volatility matrix a is typically used only to model the 
covariance matrix C for the log returns of different stocks. However, a covariance matrix 
C does not uniquely define a a such that C = era' . As we will see below, given a covariance 
matrix, the present model allows the investor to impose views on the market by selecting 
an appropriate volatility matrix that leads to expected returns of the stocks that are 
consistent with such views. Hence expected returns and investor views can be handled 
in a way that is less sensitive to statistical estimates and guesses of numerical quantities 
than existing methods. 

To explain these ideas, we first describe a basic example. It is then shown that all 
volatility matrices that imply the same C can be written as the Cholcsky decomposition 
of C multiplied by an orthogonal matrix. 

Example 2.1. Consider the two stocks S\ and S2 with covariance matrix of the log 
returns 



Assume that the companies are of about the same size and importance. It is natural to 
make the interpretation that each stock has a unique Brownian motion associated with 
it that represents the uncertainty primarily due to that stock. A consequence of this 
perception is that the volatility matrix for these two stocks should be symmetric so that 
5*i depends as much on Brownian motion B\ as S2 depends on Brownian motion £?2- 
Symmetry can be attained by taking the matrix square root of the covariance matrix. 
Given C, the matrix square root volatility matrix 




(4) 




Equation (3) gives that the continuously compounded expected returns are 



t4=r + 2.4250 - r) 
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and 



(4 = r + 2.668(^t-r). 



These quantities are similar, which is to be expected since we assume symmetric depen- 
dence between the stocks. However, if Si were considerably larger than 5 2 , it would be 
reasonable to assume that Si was independent of Brownian motion B2 , as B2 represents 
the risk primarily due to stock 6*2. Hence the volatility matrix would be the diagonal 
Cholesky decomposition, which for C is 



f4 = r + 3(fi-r). 

The difference between \i\ and [i\ is now larger, reflecting the change of volatility matrix. 

We continue with a technical lemma that helps us construct volatility matrices that 
are more advanced than the ones in the example above. 

Lemma 2.1. Assume that a covariance matrix C is given. Any volatility matrix V that 
satisfies C = W can be written as the Cholesky decomposition of C multiplied by an 
orthogonal matrix. 

Proof. We know by QR factorization that V can be written as V = LQ, where L is lower 
triangular and Q is orthogonal. It follows that C = VV = LQQ'L' = LL', regardless of 
orthogonal Q. But since L is lower triangular, it must be equal to the unique Cholesky 
decomposition of C. □ 

The implication of Lemma 2.1 is that the present model includes substantial flexibility 
in the expected returns of the stocks because, for every orthogonal matrix we use to 
rotate the Cholesky decomposition of C, we get a volatility matrix that in turn leads to 
a different set of continuously compounded rates of return for the stocks. 

One approach to finding appropriate volatility matrices is the use of orthogonal matri- 
ces to rotate the Cholesky decomposition so that the rotated matrix fits a specified target 
matrix as closely as possible in a least-square sense. The target matrix can be chosen 
arbitrarily, without considering the covariance matrix. This means that the investor can 
design freely a target matrix that, when considered as a volatility matrix, has row sums 
that lead to continuously compounded rates of return that reflect the investor's market 
views. 




In this case, 



H\ = r + 20u 



r) 



and 
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3. Markowitz' problem in continuous time 

We derive in this section an explicit solution to Markowitz' problem in continuous time, 
given our market model. The optimal strategy implicitly is to keep 1/n of the wealth 
invested in each Brownian motion, and not in each stock as in the classical 1 jn strategy. 

We define the admissible strategies A to be the set of all R™-valued stochastic processes 
that are uniformly bounded and progressively measurable in Tt ■ Note that this definition 
allows negative positions in the stocks. 

The self-financing wealth process W 71 is defined as 

'•■wrw^,., , /"< i -Er.i"<(»))w^i) dRWi 



^ (() = m + W'«£) dS , (s)+ / 

. . Jo ^i(S) Jo 



R(s) 



for all t <E [0,T], where 7Ti(i)W^ (t) / Si(t) is the number of shares of stock i held at time 
t. See [4] for a motivating discussion. We sometimes write W for W* when there is no 
risk for confusion. This gives the wealth dynamics 

n / n \ 

dW(t) = n (t)W(t) [(/" ~r)dt + dBj (t)] + W{t)r dt 

i=i \j=i / 

= W(t) I Pi (*) 0" - r) dt + Pi (*) dflj (t) + r J 
\i=i 3=1 / 

for the processes ft := ^i^ij- 

The assumption that 7r is uniformly bounded implies that the equation for W can be 
written as 



W(t) = w exp ( ^ y (s) (a* - r) - -p) (s) ds + j Pj (s) dBj (s 



rt). (5) 



We now state and prove the main result of the paper. 
Theorem 3.1. The unique solution to the continuous-time Markowitz problem 

min{Var(W w (T))}, 

TreA 

E[W W (T)] > wexp(AT), 
where A > r, is given by the strategy tt* that solves 

1 A-r 

TTlOXl H h7T n <7 n l = , 

n fi — r 



(0) 
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7T101 „ + h 7T„<7„„ = 



n (j, — r 



for all t £ [0, T] . Here w is the initial wealth, A is the continuously compounded required 
rate of return, and T is a deterministic time. Moreover, the equation for the optimal 
wealth process W* becomes 



1 f\-r\ 2 \ 1 X-r 



vr (tj^cxp^A--^— -) y + ^—Bit) 

for a Brownian motion B, where "=d" denotes equality in distribution. 

Proof. We assume that A > r such that we need to invest in some risky asset to obtain an 
expected yield larger than u>exp(AT). Further, it is a necessary condition for an optimal 
strategy that Y^j=iPj(t) — f° r all (t,u>) £ [0, T] x il. The reason is that whenever a 
strategy it violates this condition, exchanging it for the strategy to put all the money in 
the risk-free asset will both increase the expected return and lower the variance. 
For the optimal strategy 7r* we must have that 

E\W*'(T)]=wexp(XT). (7) 

To see this, consider a strategy it with E[M /7r (T)] > wexp(AT). We know that 

Vax(W n (T)) 

>E[Vax(W*(T)\p)] 



w 2 E 



cxp ^2 ^(p -r)J2 J%&) dt + rT ^j j ( CX P (E I 



pj(t)dt 



3- 

Hence, the strategy air, for any a £ (0, 1) such that 

E[M /Q7r (T)] > wexp(AT) (8) 

holds, has lower variance than n by the definition of the pj. Henceforth, we consider only 
the strategies ir £ A that satisfy equation (7). 

Set I p := -^f X)j=i Jo Pj(t) dt- We can apply Jensen's inequality to see that 



cx p(E / 

\i=i Jo 



p 2 (t)dt) >cxp(nT/ p 2 ), 



for all u) £ fi, so 

Var(M/ 7r (T)) > w 2 E[cxp(2((/i - r)nl p + r)T)(cxp(nT/ 2 ) - i)]. 

Now, set 

X = cxp((/i — r)nTl p ), 
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and note that X> 1 for all (t,ui) £ [0,T] x SI by the definition of I p . We see that 
exp(2((/i- r)nl p + r)T)(exp(nT/ 2 ) - 1) 



= exp(2rT)Af 2 (exp [nT (log (oxp(|^/ P j ) 

= cxp(2rT)X 2 (exp(nr(log(X 1 /«' i - r ) nT ))) 2 ) - 1) 

(log(AQ) 2 
(/i — r) 2 nT / 



: exp(2,T)X 2 (exp( 7 -^)-l 



But straightforward calculations show that the function 

exp(2rT)X 2 ( cxp 



/ QogPO) 2 
\ (/i — r) 2 nT 

is strictly convex in X on [l,oo), and 

E[X] = E[exp(0 - r)nTI p )} = cxp((A - r)T), 
since tt satisfies equation (7). Hence we can use Jensen's inequality to see that 

X&v(W*{T)) > w 2 cxp(2AT) (cxp ( (- — — nT 
with equality only for the strategy it* for which 

P m = ---=Pn(t) = - — , (9) 

n /i — r 

for all t £ [0,T]. Note that E[W** (T)) = wcxp(XT). Hence, for sufficiently high bounds 
on the admissible strategies, the strategy ir* that solves equation (6) for all t £ [0, T] 
minimizes the variance of the terminal wealth W(T) subject to the growth constraint in 
equation (8). Plugging in the process p* in equation (5) gives that 

A-r\ 2 t 1 X-r 



W"\t)=wexp[Xt-( — + J2 B ^) 



1 /A-r\ 2 \ 1 X-r 



= dW cx P ^A--^— j J t+ —— B (t) 
for a Brownian motion B. □ 



One intuitive approach to applying 7r* is to plug in a choice of A and estimates of 
the model parameters a, fi and r into equation (6). However, there exists an alternative 
approach. 
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n=25 
n=50 
n=100 
n=200 



0,8 1 1.2 1.4 1.6 1 2 3 4 

Figure 1. Left: Distributions of optimal wealth (1) for different n with parameters w = 1, 
A = 0.1, /i = 0.2 and r = 0.03. Right: Distributions of optimal wealth W w (1) for different n with 
parameters w = 1, A = 0.3, fi = 0.2 and r = 0.03. 

Assume that the investor has estimated a volatility matrix and chosen the total fraction 
of wealth to invest in risky assets, £™ =1 7Ti. These choices, together with equation (6), are 
sufficient to determine the strategy ir* , and the ratio (A — r)/(/i — r). This is because 
we get an additional equation by choosing a value for £™ =1 7Ti that allows us to uniquely 
determine the unknown (A — r)/(fi — r). In other words, there is no need to estimate 
either of the parameters [i and r, or to choose A, to apply the optimal strategy ir* . This 
is a considerable advantage with the present model from an applied perspective. For 
example, it is very common that institutional investors and equity funds are required to 
keep a constant market exposure, which typically is to be fully invested. Hence, for such 
investors the optimal strategy n* is uniquely determined once the volatility matrix has 
been estimated. 

The effect on the wealth process from increasing the number of stocks n is il- 

lustrated in Figure 1. The figure shows that the higher the expected return the investor 
requires, the more the investor will have to risk. Nonetheless, the risk will decrease as 
the number of stocks n increases. Note that W* is strictly positive with probability 1, 
so the investor does not risk bankruptcy. 

4. An empirical study of 7r* 

We investigate the long-term efficiency of it* by analyzing a large data set of daily prices. 
Throughout this section, the continuously compounded interest rate r = 0.03. Further, 
we assume that the investor is fully invested in the stock market at all times. Hence 
there is no need to determine the parameters fj, and A once a volatility matrix has been 
estimated, as discussed in the previous section. 

The data set, which can be downloaded from Professor Ken French's homepage, con- 
sists of the daily closing prices of 47 value-weighted industry sector portfolios. The health 
sector is removed from the original 48 sectors due to missing data. Each stock traded at 
NYSE, AMEX and NASDAQ is assigned to one, and only one, of the sector portfolios. 
In this analysis, the industry sector portfolios are treated as stocks. The data is from 
the time period 1963-07-01 to 2008-07-31. We use a "rolling-sample" approach. Every 
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Industry sector portfolios 

Figure 2. Solid and dashed lines represent the Sharpe ratio for tv* and the classical 1 jn strategy, 
respectively. The stars are the Sharpe ratios for the individual industry portfolios. 



twentieth trading day, a covariance matrix for the daily returns of the industry portfolios 
is estimated from an estimation window of the previous five years of data. The week 
of the Black Monday of 1987 is removed in the estimation of the covariance matrices, 
but it is not removed as return data. The estimated parameters are used to calculate 
the optimal strategy tt* , which in turn allows us to determine the daily returns for the 
next month. The portfolio weights are adjusted daily. This approach yields a series of 
out- of -sample returns. It is reasonable to assume that the industry portfolios arc approx- 
imately equally important to each other. Also, the investor has no preferences regarding 
any assets. Hence, we apply the matrix square root to each covariance matrix to get a 
volatility matrix, see Example 2.1. We measure performance by the out-of-sample Sharpe 
ratio, defined as the sample mean of the daily out-of-sample excess returns divided by 
their standard deviation. For this data set, the strategy tt* outperformed all but one of 
the underlying market assets in terms of Sharpe ratios; see Figure 2. Further, 7r* obtained 
15% more wealth than the classical 1/n strategy, and with 16% lower volatility. Conse- 
quently, Memmel's corrected Jobson and Korkie test of the hypothesis that the classical 
1/n strategy gives a higher Sharpe ratio than tt* had a p value smaller than 10 -4 . This 
result seemed robust in the sense that the same test with different estimation windows 
for the covariance matrix, and different frequencies with which the covariance matrix and 
the portfolio weights were updated, still yielded very small p values. Figure 3 shows the 
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Figure 3. The optimal strategy n* for the 47 industry portfolios. 



evolution of the estimated strategics 7r* , which are quite stable over time. The industry 
portfolio with the largest average fraction of wealth invested in it is the Paper index. 

We have applied the approach described above to several different data sets, although 
none that span as many years as the data set in this example. In our experience, the 
strategy it* with no investor preferences beats the classical 1/n strategy in terms of 
terminal wealth in the majority of cases, but by no means always. However, 7r* usually 
obtains lower volatility for the associated wealth process than the classical 1/n strategy. 
This typically results in a higher Sharpe ratio for , in particular when the number 
of assets is large. If various preferences or market views were used to choose a specific 
volatility matrix, naturally the result depended on how good this information turned out 
to be. 
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